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We study the vortex lines that are a feature of many random or disordered three-dimensional 
systems. These show universal statistical properties on long length scales, and geometrical phase 
transitions analogous to percolation transitions but in distinct universality classes. The field theories 
for these problems have not previously been identified, so that while many numerical studies have 
been performed, a framework for interpreting the results has been lacking. We provide such a 
framework with mappings to simple supersymmetric models. Our main focus is on vortices in short- 
range correlated complex fields, which show a geometrical phase transition that we argue is described 
by the CP k ^ k model (essentially the CP 71-1 model in the replica limit n — > 1). This can be seen 
by mapping a lattice version of the problem to a lattice gauge theory. A related field theory with a 
noncompact gauge field, the c NCCP k ^ k model', is a supersymmetric extension of the standard dual 
theory for the XY transition, and we show that XY duality gives another way to understand the 
appearance of field theories of this type. The supersymmetric descriptions yield results relevant, for 
example, to vortices in the XY model and in superfluids, to optical vortices, and to certain models 
of cosmic strings. A distinct but related field theory, the RP 2l \ 21 model (or the i?P n_1 model in the 
limit n — ► 1) describes the unoriented vortices which occur for instance in nematic liquid crystals. 
Finally, we show that in two dimensions, a lattice gauge theory analogous to that discussed in three 
dimensions gives a simple way to see the known relation between two-dimensional percolation and 
the CP k ^ k sigma model with a term. 



I. INTRODUCTION 

Ensembles of random geometric objects lead to some of 
the most subtle and interesting questions in critical phe- 
nomena. A geometric viewpoint is often useful in ana- 
lyzing phase transitions: for example, one can view an 
Ising transition as a proliferation of domain walls, or an 
XY transition as a proliferation of vortices, leading to im- 
portant dual descriptions [IH5]. But in addition to this, 
some of the most intruiguing issues concern the critical 
behaviour of the geometrical objects themselves. In this 
paper we consider the universal statistical properties of 
vortex lines in random media in three dimensions, a sub- 
ject which has received considerable numerical attention 
[6HTU], but which has so far lacked a field theoretic for- 
mulation. 

Geometrical critical behaviour can occur even in the 
absence of a conventional thermodynamic phase transi- 
tion. The best-known example of this is percolation, in 
which nonlocal geometrical correlation functions - such 
as the probability that two sites lie in the same cluster 
- yield universal exponents, despite the fact that cor- 
relation functions of local observables remain trivial at 
the critical point [11]. In fact it is precisely because of 
the short-range correlations of the local degrees of free- 
dom that the percolation universality class is ubiquitous: 
generic random media are short-range correlated. 

Two-dimensional percolation can be thought of as a 
fluctuating soup of loops [12 . On the lattice, these are 
the cluster boundaries; in the continuum, they are the 
zero lines of a random real height field v(x). When v(x) 
is biased towards either positive or negative values, the 
loops have a finite typical size, but these cases are sepa- 
rated by a critical point at which this lengthscale diverges 



and the loop ensemble becomes conformally invariant. 

Vortex lines in three dimensions constitute an analo- 
gous problem. The most generic version of it concerns 
the statistics of the zero lines of a short-range correlated 
complex field w(x). Again, these vortices form a fluc- 
tuating soup of oriented loops, with the orientation of a 
loop given by the sense in which the phase of w(x) rotates 
around the vortex. Although the correlations of w(x) are 
trivial by construction, there is an interesting geometri- 
cal phase transition as a function of bias - the ratio of 
the average value of w(x) to the width of its distribution 
[6j [10]. As will be explained in detail in the following, 
there is an extended phase at small or zero bias [EHS], 
in which infinite vortex lines proliferate (and are found 
numerically to be Brownian), while at large bias there is 
a localized phase in which long vortices are exponentially 
suppressed. In between lies a critical point, where loops 
have a nontrivial fractal structure (and a fractal dimen- 
sion close to 5/2). Similar transitions can also occur in 
the absence of applied bias, for example in the vicinity of 
the critical point of the XY model, or for line defects in 
other order parameters: in particular, one can consider 
the unoriented (Z2) vortices in nematics pTU| H3]. 

Universal aspects of vortex statistics have been exam- 
ined numerically in a wide range of contexts, including 
cosmic strings, beginning with the simulations of Vachas- 
pati and Vilenkin [7, 8| [TQ | [T3HT5] , line defects in random 
light fields, known as 'optical vortices' [9] [16] . vortices 
in the XY [HI [18] and Abelian Higgs models [EH22], 
in turbulent superfluids [23], and in a Potts-like model 
introduced to describe polymers [6]. In addition, simi- 
lar scaling behaviour has been observed in simulations of 
lattice models of completely-packed loops [24] [25] . 

Despite this previous work on vortices, numerous ques- 
tions remain unanswered, because of the lack of a contin- 
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uum formulation in which the relevant correlation func- 
tions can be expressed. In particular, it has not been 
clear how to classify these problems, or what the crite- 
ria are for two instances to be in the same universality 
class. It has been uncertain, for example, whether criti- 
cal exponents for the geometrical vortex transition coin- 
cide with those for conventional three-dimensional (3D) 
percolation, or whether universal behaviour for oriented 
and unoriented vortices coincide (even in the extended 
phase), or how the geometrical transition for vortices in 
the XY model relates to the thermodynamic one. Even 
basic facts such as the appearance of Brownian statis- 
tics in the extended phase have not been derived, though 
previously this has been explained heuristically [8] by 
analogy with polymers in the melt, long known to be 
Brownian [26 . 

Our effort to address these questions starts from a set 
of new duality transformations. We begin by showing 
that the partition function for tricolour percolation, a 
Potts-like lattice model used for the numerical investi- 
gation of vortex statistics [6] |HJ [10] , can be mapped ex- 
actly to that of a supersymmetric lattice gauge theory 
via a high temperature (graphical) expansion of the lat- 
ter. By considering the continuum limit of this lattice 
gauge theory, and, separately, by discussing mappings 
from vortex problems that are formulated directly in the 
continuum, we arrive at field theories describing the ge- 
ometrical properties of vortices. 

For oriented vortices, these field theories involve a com- 
plex supervector = (z 1 , . . . z /c+1 , x 1 , . . . x fe ), with k + 1 
bosonic and k fermionic components, coupled to a U(l) 
gauge field A. (Super symmetry [27, 28 is a standard tool 
in problems involving loops or polymers, allowing repre- 
sentation of geometrical correlation functions.) In order 
to describe the transition for vortices in short-range cor- 
related complex fields w(x), the gauge field A must be 
taken to be compact, leading to a model in the univer- 
sality class of the CP k \ k sigma model [29]. Heuristically, 
this can be thought of as the replica limit n — » 1 of the 
(jpn-i s ig ma model [30] [31], which at n = 2 is equiv- 
alent to the classical Heisenberg (or 0(3)) model. The 
ordering transition in the sigma model describes the ge- 
ometrical transition in the vortex ensemble. CP k \ k field 
theories have been related to lattice models of completely 
packed loops in both 2D - by Read and Saleur [29] and 
Candu, Jacobsen, Read and Saleur [32 - and 3D [33] [34]. 
(Strictly speaking we are abusing terminology, since there 
is a separate CP k \ k model for each integer k > 0, the 
model with k = being trivial. The correspondence with 
vortices holds for arbitrary /c, but the range of geomet- 
rical correlation functions which can be expressed in the 
theory grows with k.) 

The CP k \ k description provides the appropriate frame- 
work for questions about vortex statistics. Immedi- 
ate consequences include the explanation for Brownian 
statistics in the extended phase, which corresponds to the 
ordered phase in the sigma model language, and qual- 
itative features of the critical point (for example, the 



fact that the fractal dimension is close to df = 5/2 is 
easily understood as a symptom of confinement in the 
gauge theory). Most importantly, it clarifies the rela- 
tions between the different universality classes, allowing 
us to classify these problems and explaining why various 
numerically studied models are in the same universality 
class and others are not. For instance, it explains the 
coincidence of exponents measured for the bias induced 
transition for vortices [6] and for the completely packed 
3D loop models [24, 35 , which also permit a mapping to 
the CP k \ k model [33]. 

In the absence of bias in the distribution for w(x), the 
gauge field in the supersymmetric theory is non-compact, 
yielding the non-compact CP k \ k or NCCP k \ k model, to 
borrow the terminology used for non-compact CP 71-1 
models. The appearance of this theory can be simply 
understood by reference to a well-known duality which 
maps the XY model to an Abelian Higgs theory with a 
complex scalar field z coupled to a gauge field [3H3 l36l - 
[39] . Under this mapping, XY vortices become worldlines 
of the particles created by (to use a 2 + 1 dimen- 
sional language). However, this standard dual descrip- 
tion does not allow us to probe the long distance be- 
haviour of vortices, because it does not allow us to write 
down the required geometric correlation functions. The 
NCCP k \ k model is a supersymmetric extension of the 
Abelian Higgs model with a richer operator content (and 
possibly a richer phase diagram). 

Unoriented vortices (e.g. in a short-range correlated 
random nematic order parameter) are described by a dif- 
ferent field theory, the RP 2l \ 21 model. This is essentially 
the replica limit of the i?P n_1 sigma model. A soft spin 
version of the CP n ~ x model in the replica limit allows 
us to clarify the relation between these two geometric 
problems, as well as the relation to 3D percolation, and 
also shows that the upper critical dimension for CP k ^ k is 
six, meaning that exponents can be calculated in a 6 — e 
expansion [40-43 . The CP 1 ' 1 model also describes An- 
derson localization in symmetry class C, which is closely 
related to classical loop models [24[ [35] [44H46] , so we ex- 
pect six to be an upper critical dimension for this problem 
too. A similar conclusion has been reached previously by 
Senthil [43], who considered the soft spin formulation of 
CP 1 ' 1 . 

As mentioned above, another important conceptual 
connection is with 2D percolation. One of the contin- 
uum descriptions of this problem is as a 

C pk\k model 

with term, as was shown via a mapping to a supersym- 
metric spin chain [29j [32] . To connect this up with our 
3D results, we give a treatment of 2D percolation as a 
lattice gauge theory (with an unconventional Boltzmann 
weight for the gauge field) analogous to our treatment of 
tricolour percolation in 3D. This leads in a simple way to 
the continuum description. A similar lattice regulariza- 
tion of CP n ~ x has been discussed by Affleck [47] . 

The structure of the paper is as follows. We begin 
in Sec. fn] with a heuristic overview of the key ideas. 
In Sec. |III| we give duality mappings for oriented vor- 
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tices in short-range-correlated random systems on the 
lattice (III A) and in the continuum (IIIB). For the lat- 



tice derivation, we use a specific model, tricolour perco- 
lation [6J [8] , an elegant realization of the vortex problem 
which may be mapped to a lattice gauge theory. Sec. 



Ill C| gives the basic consequences of the field theory for 
vortex statistics. 

Sec. |IV| considers variations of the vortex problem. 
First we briefly describe the soft spin version of the 
(jpn-i moc [ e i i n the limit n —> 1 and the 6 — e expan- 



sion (Sec. IV A). We then consider a perturbation lead- 



ing to crossover to conventional percolation (Sec. IV B). 
Here we also discuss the issue of vortex intersections 
which is important in many realistic vortex problems. 
Sec. |IV C| considers unoriented vortices, occurring when 
the fundamental group of the order parameter manifold 
is Z 2 . Sec. |IV D| considers the stability of the universal 
behaviour of vortices in a short-range correlated w(x) to 
the introduction of long-range correlations. We then con- 
sider vortices near the (thermodynamic) critical point of 
the XY model, their description in terms of the NCCP k \ k 
model (Sec. IV E), and the phase structure of this field 
theory. This concludes our discussion of vortices. In 
Sec. [V] we consider loop models in 2D and their relation 
to lattice gauge theories with terms. We close with a 



short summary in Sec. VI 



II. HEURISTIC CONSIDERATIONS 

We begin with a handwaving argument as to the form 
of the field theory for vortices in a short-range-correlated 
random complex function w(x). This is intended to pro- 
vide an overview and to supplement the formal deriva- 
tions of Sec. Mil 

To reiterate, we expect the following phase structure. 
When the probability distribution for w(x) is U(l) sym- 
metric, or only weakly biased, vortices have a fractal di- 
mension of two and a finite fraction of the vortex density 
is in infinite vortex lines. This is the extended phase. A 
nonzero (translationally invariant) mean value suppresses 
vortices; when the mean is sufficiently large compared 
with the typical size of fluctuations, we enter the local- 
ized phase, and infinite vortices disappear [6j[l5]. At the 
continuous transition between these phases, vortices are 
random fractals. 

In general, a simple way to relate ensembles of loops to 
field theory is to view the loops as worldlines of quantum 
particles - in this case in 2+1 dimensions - and to ask 
what sort of particles these must be and what kind of 
interactions they must have. The partition function for 
the field theory is a sum over histories of the quantum 
system: expressed in the right basis, this is a weighted 
sum over worldline configurations. The Lagrangian must 
be chosen so that this sum reproduces the weighted sum 
over loop configurations we started with. 

Of course, to be useful the correspondence must go be- 
yond the partition function: we should be able to repre- 



sent geometrical correlation functions - such as the prob- 
ability that two points x and y lie on the same loop, 
denoted G(x — y) - as correlation functions of local op- 
erators in the field theory. The field theory we seek must 
allow us to do this despite the fact that G(x — y) is nonlo- 
cal in the original description, in the sense that we cannot 
necessarily tell, by looking only in the vicinity of x and 

whether the two points lie on the same loop. 

Since vortices have an orientation determined by the 
sense of rotation of argw(x), it would be natural to as- 
sume that we will need a single species of charged bo- 
son - with vortices that go forward and backward in the 
imaginary time direction corresponding to worldlines of 
particles and antiparticles respectively. As usual, these 
bosons would be represented by a single complex scalar 
field z(x). However, this construction would not allow us 
to write correlators like G(x — y). To be able to do that, 
we must give the worldlines (and particles) an additional 
colour label a = 1, n, so that z expands to a vector z a . 
This yields operators such as z 1 (x) z 2 * (x) , which, when 
inserted into a correlation function, forces a meeting at 
x between an outgoing worldline strand of colour 2 and 
an incoming strand of colour 1. Therefore, neglecting 
for now the possibility of strands escaping to infinity, the 
two point function (z 1 z 2 * (x) z 2 z 1 * (y)) forces x and y to 
be joined by a loop, one arm of which is of colour 1, and 
the other of colour 2 (see Fig. [T]). This is just what we 
need to construct G(x — y). By contrast, if we had only 
a single species of particle at our disposal, we could not 
separate out configurations with a loop joining x and y 
from those with x and y lying on separate loops. 

But while the colours solve this problem, they intro- 
duce another one. The sum over colour indices for each 
loop modifies the weight for a given configuration of loops 
by the unwanted factor n no ' 1o °p s : we are no longer de- 
scribing the problem of vortices, but a completely differ- 
ent problem where loops have a fugacity equal to n. 

One way to get around this is to take a replica-like limit 
n — » 1 at the end of any calculation [48 . A more concrete 
alternative [27-29 is to exchange z for a supervector 



1> = (z 1 ,... 



x\-x k ) 



(1) 



with k + 1 bosonic and k fermionic components. 
Fermionic loops come with a minus sign. Thus, summing 
over the possible colours for a loop yields a fugacity of 






FIG. 1. (Color online) Left: The two-point function 
(z 1 z 2 * (x)z 2 z 1 * (y)) forces x and y to be joined by a loop. 
With only a single colour index, we could not separate such 
configurations from those with the topology shown on the 
right. 
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k + l — = 1 as desired, independently of k. 

What should the Lagrangian for i/? be? It must be in- 
variant under supersymmetry transformations preserving 
i/j^i/). At a technical level, this ensures that the partition 
function is independent of the number of fermions k. In 
terms of loops, the invariance is a consequence of the fact 
that the colour index we have introduced is just a label: 
the weight of a configuration does not change if we alter 
the colour of one of the loops (except by minus signs as- 
sociated with fermionic colours). It is tempting also to 
assume that since w(x) is short-range correlated, interac- 
tions between vortices will be short-range, and therefore 
that we should write down a field theory with local in- 
teractions. The natural candidate is 

incorrect = | V</f + /i|^| 2 + A|^| 4 . (2) 

However, this is not what we want. The phase transition 
(as a function of /i) at which xj) condenses does indeed 
represent a proliferation of rj) worldlines. But this phase 
transition is thermodynamically nontrivial - by setting 
k = we see it is in the XY universality class - and 
so cannot represent the vortex phase transition we are 
discussing, where all local degrees of freedom are short- 
range correlated. 

As a hint towards a better candidate, notice that cor- 
relation functions of the form 

(Mx)] N *[r(x)] N °« (...)), (3) 

where (...) contains insertions away from x (and we omit 
colour indices) have no meaning in the vortex problem 
unless 7V in = Af out : for obvious topological reasons, the 
number of vortex strands entering the vicinity of x must 
equal the number leaving. We are thus free to set these 
correlation functions to zero, which can be done by in- 
troducing a gauge symmetry: 

£=|(V-a)VH 2 + ^|V| 2 + A|V>| 4 + -.- (4) 

The significance of this change depends on what La- 
grangian for the gauge field is hidden in the ellipsis in 
Q. In the trivial limit of zero gauge coupling, when 
non-gauge fluctuations of A are completely suppressed, 
the thermodynamics of Q will be unchanged and we 
will not have made any progress. On the other hand, a 
nontrivial action for A will mediate interactions between 
particles which depend on their charge. This is what 
we want - it reflects the fact that the weight of a given 
configuration of vortices depends on their relative orien- 
tations, as is clear from thinking about the winding of 
3igw{x). 

Our previous assumption of local interactions for -0 
was in fact incorrect. The weight for a vortex configu- 
ration C comes from integrating over all compatible con- 
figurations of w, and this yields a long-range interaction 
between vortices. The role of the integral over A is to 
simulate this part of the statistical weight. This is the 
fundamental idea of XY duality G3EJ l36ti39| [49] - in the 
context of the present discussion, see especially [4| [38]. 



The coupling of loops to a gauge field has also been ex- 
ploited to derive exact results for self-avoiding loops in 
two dimensions [50j [ST] . 

Specifically, in the presence of the gauge field each con- 
figuration C of oriented loops is weighted by the expec- 
tation value 

W{C) = (expi J AAx^j , (5) 

of a Wilson loop calculated using the functional integral 
over A only, and it is the desired form of this weight 
which determines the required action for A. It is par- 
ticularly easy to see what this should be in the limit of 
large bias. Suppose that (w) is large and positive, and 
consider W(C) for a single large loop. Then argiu(x) will 
be close to zero for most x, but the vortex must be the 
boundary of a sheet where aigw(x) = tt. This sheet will 
cost an 'energy' scaling as the minimal area A enclosed 
by C, so that log W(C) ~ -A. 

Such an area law for the Wilson loop implies that A 
should be a compact gauge field [38 , 52, 53 : we must in- 
clude singular Dirac monopole configurations in the func- 
tional integral over A. While we used large bias for this 
discussion, monopoles are in fact present at any nonzero 
bias. (Since they are irrelevant in the extended phase, 
the zero-bias point is not distinguished from the rest of 
the extended phase as far as the universal statistics of 
vortices are concerned.) 

The field theory we have arrived at, with a supervec- 
tor coupled to a compact gauge field A as in Eq. Q, will 
appear naturally in our derivations. However it is often 
more convenient to use alternative formulations which 
describe the same universal physics. We now discuss two 
such formulations. We temporarily simplify the discus- 
sion by exchanging tj) for n-component bosonic vector z, 
with the replica limit n — » 1 in mind. 

For the first, we switch to a sigma model, with the 
constraint z^ z = 1 rather than a potential for z, and the 
Lagrangian: 

C CP ^ = \ (|V*| 2 - |* f Vz| 2 ) , ztz = 1. (6) 

The peculiar-looking kinetic term is designed to be in- 
sensitive to the phase of z, so that we retain the gauge 
symmetry z — >• e^z. Because of this gauge symmetry, 
the local degrees of freedom live on CP n_1 , the mani- 
fold of n-component complex unit vectors z subject to 
the identification z ~ e^z [3(3 [31]. The sigma model 
describes the ordering transition of 'spins' living on this 
space. At long distances it is equivalent to the compact 
gauge theory, in which the effect of the gauge field is to 
confine the charged particles z into neutral composite 
particles, or alternately to render the overall phase of z 
a redundant degree of freedom. 

The order parameter for the transition in this sigma 
model is a gauge- invariant traceless n x n matrix, 

Q = zz ] - (z ] z)/n, trQ = 0. (7) 
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In the limit n —> 1, or in the super symmetric formulation 
described below, the transition is that between localized 
and extended vortex phases, and the extended phase is 
that in which (Q) ^ 0. The phase diagram is shown in 
Fig.0 

When n = 2, the order-disorder transition is simply 
that of the 0(3) model, to which the CP 1 model is 
equivalent. This equivalence is seen by using the Pauli 
matrices to form a three-component Heisenberg spin, 
S % = tr (j % Q = z^o % z\ when written in terms of 5, the 
sigma model Lagrangian (|6| becomes that of the 0(3) 
sigma model. The fact that n = 1 (the 'CP ' model) 
must be regarded as a limiting case is seen in the fact 
that Q = in this limit. 



Bias 

Extended Critical 



Localized 



<Q>*0 



<Q>=0 



FIG. 2. (Color online) Phase diagram for vortices in short- 
range correlated w(x). 

A familiar idea in the context of the 0(3) model is that 
a soft spin formulation in which S 2 is unconstrained is 
often more useful than a sigma model in which S 2 = 1. 
The second formulation we will need is an analogous soft 
spin theory for Q, 



£ so ft spin = tr ( VQ) 2 + fi tr Q 2 + v tr Q 3 



(the cubic term vanishes in the special case n = 2). This 
formulation, in the n — » 1 limit, will be useful to us in 



Sec. [TV] in considering the effect of perturbations to the 
vortex problem. An immediate inference from it is that 
the upper critical dimension for this universality class 
(and for the related problem of the Anderson transition 
in class C [43]) is six. Exponents for vortices can thus 
be computed in a 6 — e expansion similar to that for 
percolation [4QH42] . 

The supersymmetric replacement for the replica limit 
of the sigma model |6| is got simply by replacing z with 

</>, 



c cpk]k = 1 (|vvf - |V f Vt/f ) , V>V = i. 



(9) 



In two dimensions, these supersymmetric sigma models 
and their connection with loop models have been studied 
in [29] [32] . The gauge- invariant operator Q becomes a 
supermatrix [54]. 



str Q = 0. 



(10) 



The transition into the ordered (extended) phase as a 
function of the coupling g is equivalent to the sponta- 
neous breaking of supersymmetry and the appearance of 
an expectation value (Q). 

This simple phase diagram can be extended in various 
ways. We can for example modify the problem so as to in- 
duce a crossover to the universality class of conventional 



percolation, or to that for nematic vortices (Sec. IV). The 
gauge theory Eq. Q also prompts us to ask about the ge- 
ometrical critical behaviour described by the version with 
a noncompact gauge field. Adopting the standard termi- 
nology for the case without fermions [55 , we will call 
this the noncompact CP k \ k model, or NCCP k \ k . Since 
compactness (the presence of Dirac monopoles) is a con- 
sequence of bias, the NCCP k ^ k model describes the case 
without an applied bias. It is a supersymmetric exten- 
sion of the well-known dual theory for the XY transition 
[3H3 136H39] , and is relevant to the statistics of vortices in 
the XY model at zero magnetic field [TT] [18] . Here we re- 
gard w(x) as the order parameter for an XY model, and 
we are no longer necessarily in the domain of short-range 
correlations. We will consider this noncompact theory in 



Sec. IV E An added complication, in comparison to the 
compact case, is that it may show two separate phase 
transitions, one thermodynamic (associated with the on- 
set of the Higgs mechanism) and one geometrical (as- 
sociated with the breaking of supersymmetry). This is 
suggested by numerical results for XY vortices [T7J [18] . 
Surprisingly, in this scenario the geometrical transition 
may again be described by a compact CP k ^ k model. 



III. VORTICES IN A SHORT-RANGE 
CORRELATED RANDOM ENVIRONMENT 



(8) We now turn to derivations of the CP k ^ k model for vor- 



tices, beginning (Sec. |III A) with a lattice model in which 



geometrical observables can be explicitly related to cor- 
relation functions of ?/?. We then present a derivation 
directly in the continuum in order to emphasize the uni- 
versality of the results and to see the r elation with more 
standard treatments of duality (Sec. IIIB|). Sec. |III C 



then gives simple consequences of the field theory de- 
scription. 



A. Tricolour percolation 

Tricolour percolation is an elegant discretization of the 
vortex problem [6j [8] which has become a standard setup 
for numerical simulations, partly because some more ob- 
vious alternatives to this model suffer from aesthetic de- 
ficiencies. For example, while it would be natural to con- 
sider random phases = argw on the sites of a cubic 
lattice, with vortex lines piercing the plaquettes, this re- 
quires an additional rule for deciding when a plaquette 
is pierced by a vortex line. Also, in such a model noth- 
ing prevents vortex lines from intersecting each other, 
requiring a further decision about how to resolve the in- 
tersections. For a generic function w(x) in the continuum 
on the other hand, zero lines intersect with probability 
zero. (Intersections will be discussed in detail in 



IV B 



) 

Tricolour percolation avoids these issues by a smart 
choice of lattice, and by discretizing not only real space 
but also the target space for the degrees of freedom. This 
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idea is familiar from 2D percolation, where a single uni- 
versality class encompasses both continuum percolation, 
which deals with zero lines of a continuous random func- 
tion v(x), and lattice site percolation, where the function 
takes only the values v = ±1, and where intersection of 
cluster boundaries is avoided by using sites on the tri- 
angular rather than the square lattice. For a random 
complex function, the minimal discretization is to let w 
run over the cube roots of unity, w = 1, e 27 ™/ 3 , e 47 ™/ 3 [7]. 
The configurations are therefore configurations in perco- 
lation with three colours [6 , and vortices (or tricords) 
are lines where all three colours meet. 

The construction described below was introduced by 
Scherrer and Frieman [8 as an improvement on the model 
of Vachaspati and Vilenkin [7 and continues to be used 
for numerics in the context of cosmic strings [10]. It 
was introduced independently and simulated extensively 
by Bradley, Debierre and Strenski [6 , who realized that 
the model displays a novel geometric phase transition 
at nonzero bias. Their interest came from a different 
direction: they were looking for an efficient means of 
generating configurations of self-avoiding polymers, and 
were inspired by earlier work on smart walks [56] . 



1. Model 

Fig. [3] depicts a tiling of space with cells of a certain shape 
(truncated octahedra). The centres of these cells lie on 
a bcc lattice, and they are the Wigner-Seitz cells for this 
lattice. The loops will live on the edges of the cells, which 
form a second lattice known as the tetrakaidekahedral 
lattice. Its key features are that each link is surrounded 
by three cells, and that its sites are four-coordinated. 

Tricolour percolation configurations are generated by 
randomly and independently colouring the cells either 
red, green or blue, with probabilities £>r, pg and p&. As- 
signing each colour to a cube root of unity 

{R,G,B} — > {1,cj,cj 2 }, u = e 27ri/3 , (11) 




FIG. 3. (Color online) The lattice for tricolour percolation. 
Left: showing the relation to the bcc lattice. Right: with the 
locations of the dual degrees of freedom (-0 on sites and U on 
links) indicated. 



gives the relation with a discrete random field w(x); the 
threefold anisotropy in its distribution is not expe cted t o 
change the universality class of the transition (Sec. |IIIB ). 
Vortices are the lines along which all three colours meet. 
The design of the lattice ensures that these curves are 
non-branching and self and mutually avoiding, and they 
have an orientation defined by the order in which the 
three colours encircle the link. 

Close to the symmetric point, £>r = pg = Pb = 1/3, 
tricords are in the extended phase [6j [10]. Moving 
sufficiently far from this point, we enter the localized 
phase. Numerical studies of the transition have been 
performed only on the line £>r ^ Pg = Pr, which in- 
tersects the critical line at the two points £>r = 0.417(1) 
and £>r = 0.255(5) [6], but the derivation below leads us 
to believe that the critical behaviour will be the same at 
all points on the critical line, as expected by the authors 
of [6]. Note that the critical line lies inside the region 
in which all three colours percolate, so the transitions do 
not coincide with percolation transitions of the colours 
(the critical probability for site percolation on the bcc 
lattice is 0.2460(3) [6]). 

2. Symmetric point 

We begin with the symmetric point, where the mapping 
to lattice gauge theory is simplest. The gauge theory we 
will need is unconventional in that its Boltzmann weight 
is not naturally written as the exponential of something 
simple. This situation is standard in the study of loop en- 
sembles: lattice magnets with Boltzmann weights of non- 
exponential form are useful because they give straightfor- 
ward graphical expansions, which can be interpreted as 
loop models [57] . For a pedagogical introduction to these 
ideas, see [58 . The motivation here is similar: a graphi- 
cal expansion of the kind standard in lattice gauge theory 
[1] will yield tricolour percolation configurations. 

Let z, j . . . label sites of the tetrakaidekahedral lattice, 
and let (ij) denote a link directed from i to j. We intro- 
duce gauge fields Uij on the links, which are unimodular 
numbers with Uij = U*^ and unit supervectors [27] [28] 

^ = (4,-,^ +1 ,x!,-,xf), 4^ = 1 (12) 




FIG. 4. (Color online) A section of tricord in tricolour perco- 
lation. 
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at sites. 'Tr' will denote an integral over all degrees of 
freedom, with the normalization Tr 1 = 1. The normal- 
ization of ip then implies that Tr ipfip*^ = 5 a P . 

Letting F denote a face (where two cells meet) and I 
a link, the partition function for the lattice gauge theory 
is 

z = ttI[Ii + IIu 1 + 1[uA 

f \ ieF leF ) 

xf](l + [/^.+c.c.). (13) 

(ij) 

The links in the product YlieF are oriented consistently 
around the face F, and changing the orientation corre- 
sponds to exchanging Yl U and Yl U* . 

To see the relation to tricolour percolation, we expand 
out the two products in Eq. ( [l3| ) - one over faces F and 
one over links (ij) - and represent the terms graphically. 
The diagram for a given term is built up as follows. For 
each face F, we must choose either 1, f|/7 or f|/7*. If 
we choose '1' we add nothing to the diagram, while if 
we choose Y\U or Yl U* we draw in the face F together 
with an orientation (equal or opposite to that of the links 
in the product respectively). We represent this orienta- 
tion by a normal vector as in Fig. [5] For each link we 
must choose either 1, in which case we draw nothing, or 
Ufji^l^j, in which case we draw in the link with an ar- 
row from j to z, or U^^tf)^ in which case we orient the 
arrow the other way. 

Only a subset of the diagrams generated this way sur- 
vive after integrating over U, namely those having an 
equal number of JJ\ and U* on each link I. For each link 
there are three possibilities (shown in Fig. |5|. Either (i) 
none of the three faces bordered by the link is included 
in the diagram, and neither is the link; (ii) two of the 
faces are included, they are consistently oriented so as 
to form part of a sheet of oriented surface, and the link 
is not included; or (iii) the link is included and so are 
all three of the faces, with orientations determined by 
that of the link and the right hand rule. We are then 
left with sheets of oriented surface which close on them- 
selves or meet along directed lines C. On a finite lattice, 
these lines are closed loops. They are self and mutually 
avoiding, and neither branch nor terminate. 

Qi} 0$ 

FIG. 5. (Color online) Possibilities for the graphical expan- 
sion at a link, as described in the text. In the mapping to tri- 
colour percolation, the shaded surfaces become domain walls 
between different colours. 



Such configurations are easily mapped to tricolour per- 
colation by regarding the sheets as boundaries between 
domains of different colour. Fixing the colour of one cell, 
we can colour all the others by the rule that the colour 
changes cyclically, 

R G -> B -> R, (14) 

upon crossing a sheet in the same direction as its normal. 
The lines C where three sheets meet are then the tricords. 
A convenient convention is to consider a finite bcc lattice 
and define the colour outside the lattice to be (say) red. 
Thus a cell on the boundary is red if its exterior faces 
(those on the boundary) are not covered by surface in 
the graphical expansion. 

Let T denote a tricolour percolation configuration with 
loops C. At this point, having integrated over U but not 
i/?, we have 

z = e I Tr n (^ 2 )....(^i) | , (is) 

where 1, . . . , L label the sites on a given loop C. To eval- 
uate the integral over the supervectors, note that each 
loop gives 

Tr (v>r mfM^V) ■ ■ ■ wtww) ■ (16) 

The formula Tr^f^f* = S a ^ requires all colour indices 
to be equal, say to a, to give a nonzero result. For the 
fields ipi and however, a reordering is necessary in 
order to use this formula: this gives a minus sign when 
a corresponds to a fermionic component of Thus the 
sum on a yields a fugacity per loop of fc + 1 — fc = 1, 
and Z becomes as a sum over equally weighted tricolour 
percolation configurations, Z = J^-l = 3 no ' cells . 

In addition, we can use i/?s to construct geometrical 
correlation functions. Let the operators Oi depend only 
on the i/?s. Then the graphical expansion goes through 
as before and we have 

(Oi . . . o N ) = 1 Tr °i • • • °* II (v4</> 2 )....w^i). 

r cer 

This formula allows us to translate correlators of i/)s into 
the probabilities of various geometrically-defined events. 
The other input required is a formula for the single site 
traces, 

Tr (ip ai ...^i/j*? 1 ...^*^) 

= ^ (ij ai ...ip^Tp*? 1 ...^ g ) G , (17) 

where (. . .) G is a Gaussian average evaluated using 
Wick's theorem and (?/; Q! ?/;*^) G = 5 a ^ . A simple example 
is the two-point function (Q^Q^ 1 ) of the gauge- invariant 
operator Qj 2 = zjz* 2 . The operator insertions restrict 
the sum over T to configurations in which i and j lie 
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on the same tricord, and so gives the probability of this 
event. Note that in order to construct this correlation 
function we require k > 1. We will return to correlation 
functions in Sec. IIII CI 

The naive continuum limit of the lattice gauge the- 
ory (13) involves a superfield coupled to a gauge field 



A. The compactness of the microscopic gauge degrees of 
freedom XJ\ translates in the continuum into the presence 
of Dirac monopoles in A [52] [53]. However, because i/j 
couples to U 3 rather than to U in (13), these monopoles 



are of three times the minimal charge allowed by the 
Dirac quantization condition. (We will take the minimal 
value allowed by the Dirac quantization condition as our 
unit of monopole charge - hence the factor of three in 
Eq. fT8l below.) 



Crudely, monopoles arise because there is a 67r ambi- 
guity in the magnetic flux Bp passing through each face 
F - here F denotes a face together with a choice of orien- 
tation - and thus also in the divergence of this quantity. 
Bp is a vector field on the links of the dual lattice, which 
pierce the faces F, and is defined modulo 6tt by 



exp 



iBi 



(18) 



leF 



The divergence (V .B)& is defined at each site of the dual 
lattice, or equivalently each tricolour percolation cell S, 
by the sum of Bp over the (outwardly oriented) faces of 
the cell. 

Recall that an ambiguity of this kind can be resolved in 
two ways, which correspond to thinking about B either 
as a vector field whose divergence is nonzero at monopole 
defects, or as a vector field with strictly zero divergence 
but with Dirac strings. If we resolve the ambiguity by re- 
stricting Bp to the domain [— 37r, 37r), it can have nonzero 
divergence at the locations of monopoles, whose charge 
p is a multiple of three. Letting B denote a site of the 
dual lattice, or equivalently a tricolour percolation cell, 



(V.B) B = 27T PB , 



pb e 3Z. 



(19) 



Alternatively we can resolve the ambiguity (call the new 
version B r ) in such a way that 



(V.B')b = 0. 



(20) 



B' differs from B by the inclusion of Dirac strings, which 
carry 6ir flux away from the monopoles along strings of 
adjacent plaquettes, and on which B' lies outside the 
region [— 37r, 3tt). 

In a continuum formulation, the functional inte- 
gral over A must include configurations with pointlike 
monopoles (of charge p G 3Z), together with singular 
Dirac strings which do not incur any weight in the ac- 
tion [53j [59] . A fugacity will also be associated with 
each monopole. Anticipating Sec. IIIB[ these charge 
three monopoles are a result of threefold anisotropy in 
the argument of our random field argu> [38] . We post- 
pone further discussion until Sec. IIIB[ since they will 



not play an important role. The symmetric point lies 
in the phase where i/? is condensed and A is massive by 
the Higgs mechanism; the appropriate description here is 
in terms of Goldstone modes associated with the break- 
ing of supersymmetry, and the monopoles are irrelevant 
[81]. However, monopoles of charge one appear on mov- 
ing away from the symmetric point, and will play an 
important role at the transition. 



3. Alternative view of the graphical expansion 

For an alternative view on the graphical expansion, and 
to connect up with the discussion in Sec. [H] we can ex- 
pand out only the second product (that over links) in the 
partition function (13). The diagrams now consist only 



of oriented links (to each of which is attached a factor 
U 3 ^^). Next, integrating over ?/j kills all configurations 
except those in which the links form closed loops. How- 
ever it does not prevent these loops from intersecting. If 
C denotes such a configuration of loops, then 

^=£Tr,n(i+n^+n^) n (m)- 

c f \ leF leF J loops c \iec / 

In the worldline picture, we see that each loop configu- 
ration is weighted by the expectation value of a Wilson 
loopWXC), 



w)=(nn^ 3 ) > ( 2i ) 

c \ceciec I v 

evaluated using the [/-dependent part of the Boltzmann 
weight: 



{■■■), 



^(..onfi+n^+n^Y (22) 

u f \ leF leF / 



The expectation value W(C) assigns the correct entropic 
weight to each tricord configuration. This is a nontrivial 
weight given by summing over the possible domain wall 
configurations compatible with a given set of tricords. In 
particular the weight is zero for loop configurations which 
are not allowable tricord configurations, including those 
with intersections of loops. 



4. Away from the symmetric point 

At first sight, extending the correspondence between tri- 
colour percolation and lattice gauge theory away from 
the symmetric point presents a problem, since the basic 
objects in the graphical expansion of the lattice gauge 
theory are the domain walls between colours rather than 
the colours themselves. But in fact this is not the case. 
To begin with we consider changing the probabilities 
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for a single cell: this allows the appearance of a Dirac 
monopole of unit charge at the centre of that cell. 

Again, we consider a system with boundary, and stip- 
ulate that the colour outside the boundary is R. The 
colour of an interior cell B is then determined by the 
signed number of domain walls we must cross (modulo 
three) in order to reach a point outside the boundary, 
starting from B. 

Pick any path V from the interior of B to the boundary. 
Now introduce a variable 77 which will run over the cube 
roots of unity, and modify the Boltzmann weight for the 
lattice gauge theory on the faces crossed by V by the 
substitution 

V l£F IczT? / V Icz T7 IczT? 



leF 



ieF 



leF 



The factor 77 goes with the term in which the plaquette's 
normal is parallel to the path V, and the factor 77* with 
the term where the normal and V are antiparallel. 

For a given term in the graphical expansion, we can 
read off the colour of B from the power of 77 mod three: 
if all the 77* s and 77s cancel, B is red, and if there is an 77* 
or an 77 left over, B is green or blue respectively. Thus to 
change the probabilities away from 1/3 for this cell, we 
make 77 a dynamical variable and include in the partition 
function a sum over 77 with the weight 



PgV+PbV*) (• 



(23) 



77=1, UJ, UJ 2 



where the ellipsis stands for the other factors in the par- 
tition function. The graphical expansion goes through 
as before, with each tricolour percolation configuration 
acquiring a power of 77 determined by the colour of B. 
Performing the sum over 77, only one of the three terms 
PR, PgV or PbV* gives a nonzero contribution, yielding 
the appropriate probability. 

We now specialize to the line £>r > Pg_=Pb (we will 
briefly discuss the general case in Sec. IIIB| ). Writing 




FIG. 6. (Color online) Away from the symmetric point, 
monopoles of unit charge appear. In the presence of 
monopoles, the Boltzmann weight for the gauge field is mod- 
ified on the plaquettes crossed by Dirac strings. 



77 = cj Pb , where p% G {0, ±1} will become the monopole 
charge in the cell, the above weight is 



PB=0,±1 



h = 



PR 



(24) 



At the symmetric point, p& is forced to be zero, and the 
modification disappears. Moving away from the symmet- 
ric point allows ps to fluctuate. 

To see the effect of this on the gauge theory, consider 
the Boltzmann weight in terms of gauge field configura- 
tions, rather than in terms of the graphical expansion. 
The Boltzmann weight for the plaquettes on V is now 
maximized not by a magnetic flux Bp = (mod 67r), 
but by a flux Bp = 2irps (mod 6tt) oriented antiparallel 
to V. Thus V is a Dirac string for a monopole of charge 
pB located at 6, and h is a fugacity for monopoles. 

We may immediately extend the above argument to a 
uniform (rather than local) change in the probabilities. 
Each cell acquires a monopole charge and we write 

Z = J2 hNpr R II i 1 + e i(BF ~ DF)/3 + c.c.) 



{PB} 



1I(i • ■ r;^>,). 



(ij) 



(25) 



Here N p is the total number of monopoles, N p = J2 B \p&\. 
Bp is the magnetic flux through the face F, defined in 
(18), and Dp is the flux through F due to the Dirac 



strings V emanating from the monopoles. Since the ge- 
ometry of these strings is arbitrary, Dp is arbitrary ex- 
cept that it must be a multiple of 2tt on every face, and 
(viewed as a vector field on the dual lattice) must have 
the right divergence, (V.D)& = —2i\p^. We fix it by any 
suitable convention, for example by setting all the Dirac 
strings to be parallel to a given axis. 

A peculiarity of the above mapping is that monopoles 
appear in two distinct ways: as a result of the microscopic 
compactness of U and as a result of a modification to the 
Boltzmann weight for the gauge field. In the continuum 
derivation (Sec. IIIB) all monopoles appear in the same 
way. 



5. Interpretation 



The lattice model (|25| describes a supersymmetric mat- 
ter field, whose worldlines are the tricords, coupled to a 
gauge field. The partition function includes a sum over 
monopole configurations, with monopoles weighted by a 
fugacity that tends to zero at the symmetric point. If we 
write a continuum Lagrangian with a gauge field, 



£ = |(V - iA)i>\ 2 + + AM 4 + kB 2 



(26) 



we must take into account that in the kinetic term B 2 for 
the gauge field we should not include contributions from 
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Dirac strings, and that the monopole fugacity is hidde n 
in the measure for the path integral [53j|59] (Sec. IIIB). 
But as mentioned in Sec. [TTJ the universal behaviour is 
equally well captured by a sigma model formulation. The 
Lagrangian for this sigma model is given in Eq. 

The microscopic form of the action (25) does not of 



course give us the values of the parameters (such as \i 
or g) in the coarse-grained action. However it is easy 
to identify the phases of tricolour percolation with the 
phases of the above theory. 

The extended phase is the ordered phase of the sigma 
model at small g, or the Higgs phase (at negative renor- 
malized /i) in the language of the gauge theory. The 
gauge symmetry (Higgs mechanism) removes one Gold- 
stone mode, but 2k bosonic and 2k fermionic Goldstone 
modes remain. They are governed by free field theory, 
and this leads to Brownian behaviour for the vortex lines 



(see Sec. IIIC for further details). Note that that while 
the symmetric point may seem special, in that charge-one 
monopoles are forbidden there, it is no different from the 
rest of the extended phase as far as universal properties 
are concerned. 

The localized phase is the disordered phase of the sigma 
model. In the gauge theory, both i/? and A are massive, 
the latter because of the proliferation of monopoles. 

Note that when k = the two phases collapse into one, 
and the Goldstone modes which distinguish the extended 
from the localized phase disappear. This reflects the fact 
that when k = we can no longer construct the geomet- 
rical correlation functions needed to probe the transition. 



B. Vortices in the continuum 

Focussing on tricolour percolation allowed us to introduce 
the dual degrees of freedom in an explicit way, but it had 
the disadvantage that the more general applicability of 
the CP k \ k model was obscured by microscopic details 
specific to that model, such as threefold anisotropy. In 
order to make the universality of the result plausible, 
we now give an alternative derivation in the continuum, 
building on standard duality arguments for the XY model 

laHa ESH39] . 

Let 9{x) denote the phase of the random complex func- 
tion, 



9{x) = argnj(x). 



(27) 



This field is compact, 0(x) = 0(x) + 27r, and ill-defined 
at the vortices. We expect universal features of their 
statistics to be independent of the precise details of its 
probability distribution, so long as correlations in Q(x) 
are short-range (just as universal quantities in percola- 
tion are robust to changes in the probability distribution 
of the microscopic variables). Thus we choose the most 
convenient distribution, taking 0(x) to be governed by 
the Lagrangian £xy of the classical XY model, 



When H = 0, this has a phase transition at a critical 
coupling K c , with short-range correlations in the disor- 
dered regime K < K c . The magnetic field H allows us 
to introduce bias, i.e. to tune the mean value of the ran- 
dom function w(x): at K < K c and sufficiently small \H\ 
vortices are in the extended phase, and we enter the lo- 
calized phase by increasing \H\. This gives a critical line 
separating extended from localized vortices. Note that 
in this formulation the core energy of a vortex is hidden 
in the short-distance regularization. 

As is well known, the phase transition of the 3D XY 
model at zero magnetic field is dual to a Higgs transition 
for a Lagrangian with a non-compact gauge field coupled 
to a single complex scalar z [3HSJ I36H39] , 

Aiuai = |(V - iA)z\ 2 + k(V x A) 2 + v\z\ 2 + \\z\\ 

where the role of the gauge field A is to encode the long- 
range interactions between vortices which come from in- 
tegrating over all 6{x) compatible with a given vortex 
configuration. The disordered phase of the XY model, 
in which vortices proliferate, corresponds to the Higgs 
phase of the gauge theory (at negative renormalized /i), 
where the field z has condensed and A is massive by the 
Higgs mechanism. The ordered phase of the XY model 
corresponds to the disordered phase of the noncompact 
dual theory in which z is massive and A massless. 

This dual theory differs from the one we require in two 
ways. Firstly, we must consider the theory at nonzero 
H. The operators e ±z0 correspond in the dual language 
to the insertion of magnetic monopoles, so the perturba- 
tion HcosQ leads to compactness of the gauge field, a 
point also stressed by Kleinert [38 . Usually this pertur- 
bation would lead simply to a massive theory. However, 
in order to capture geometrical correlation functions, we 
must extend the theory, using supersymmetry to intro- 
duce extra degrees of freedom. This allows us to obtain 
a nontrivial field theory (the CP k \ k model) at finite H. 

Following the general approach described in [39], the 
functional integral over 0(x) can be split into a 'sum' over 
vortex configurations C, and for each C, an integral over 
compatible field configurations. To do this, we associate 
a representative field configuration <p(x) with each vor- 
tex configuration (by any suitable convention), and write 
6{x) in terms of <p(x) and a vortex-free field Q(x): 



o(x) = </>(x) + e(x). 



(29) 



6{x) should be treated as a single-valued field in 



£ XY = ^(V6>) 2 -HcosO. 



(28) 



[—00,00): while Eq. (29) defines 6{x) only modulo 27r, 
this ambiguity can be removed using the line integral of 
V[0 -</>]. 

The path integral over now becomes J c J V6, where 
J c is a shorthand for the sum over vortex configurations. 
To make this precise we would have to fix a regulariza- 
tion of the theory. This will not be necessary for our 
present purposes; however it is important to note that 
J c includes a local weight determined by physics on the 
scale of the UV cutoff. This will include an action cost 
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per unit length of vortex, and local interactions between 
vortex strands. Note that one way to think of the sigma 
model action (28) is as an extreme case in which fluctu- 
ations in the amplitude of w ~ e %e are restricted to to 
vortex tubes with a thickness on the scale of the short- 
distance cutoff. 



Inserting the expression ( 29 ) for into the action, and 
making a Hubbard-Stratonovich transformation with a 
three-component vector field £ [39], one obtains the La- 
grangian 



XY 



■«f.(V0 + V0) -HcosO. (30) 



When H = 0, the integral over forces £ to be 
divergence- free. A nonzero H relaxes this constraint, al- 
lowing point-like monopole defects where V.£ ^ 0. This 
can be seen by a standard trick [53] - expanding the 
Boltzmann weight in H gives 

exp(i//d3,cos^))=^exp( i /d3,^ 

(31) 

where p{x) is a density for a variable number N of 
monopoles, of strengths pi = ±1, at locations r^, 

N 
i=l 

and J p is a shorthand for an integral over monopole con- 
figurations, weighted by a fugacity H/2 per monopole 
and with a factor of 1/N\ for indistinguishability: 

/ denotes £ ^ ft f E I ^) . (32) 

J P N=0 i=l \ P i=±l J J 

Integrating over the XY partition function is 
Z XY = [ / / V£5(V4-P) 

JC Jo J 



ICJp 

x exp 



i 2 + i£.V<£ + ip<t> 



(33) 



Since the point-like charges in p are quantized, the 2tt 
ambiguity in <fi does not affect the Boltzmann weight. 

We resolve £ into a divergence-free part V x A and a 
collection of Dirac strings D which carry the flux away 
from the monopoles along singular flux lines, by writing 

^=^(Vxi-D), with V.D = -2irp. (34) 
2tt 

(The geometry of the Dirac strings is chosen by some 
convention and does not fluctuate.) 

As a consequence of (34), the functional integral over 



A includes Dirac monopole configurations with singular 
strings D. The form of the kinetic term for the gauge 
field, involving £ 2 rather than (V x A) 2 , means that these 



strings do not cost any energy. Actions of this form for 
compact gauge fields are discussed in [53j [59] . 

We now simplify the action in Eq. (33 ). Since V(j) is not 



the derivative of a continuous, single- valued function (in 
terms of such a derivative, V0 = — ze _ *^Ve^), the last 
two terms do not cancel on integrating by parts. Instead 
they yield: 



;|d 3 x(c. 



+ P 4>) 



~2tt 



- I d 3 iA(VxV^ + ! f D.d 2 S. (35) 

7T J J surfaces of 

discontinuity 



The second integral counts intersections of the Dirac 
strings with the surfaces on which <j) jumps by 2tt. Since 
this is an integer multiplied by 27r, the term does not 
contribute to the Boltzmann weight. However, the first 
term provides a coupling of the gauge field to a singular 
current 



J -^V x V0 

Z7T 



(36) 



which runs along the vortices. 

At this point we have a theory of fluctuating vortex 
lines - flux lines of J - coupled to a compact gauge field, 
with partition function 



Note that because of the quantization of the monopole 
charge and of the vortex 'flux', the Dirac strings do not 
couple to J. The dependence on H is in the integral over 
monopole configurations (32) - when H = 0, monopoles 



are absent and A is noncompact. 

The integral over vortex configurations, so far denoted 
only symbolically by J c , is over configurations of ori- 
ented loops, with local interactions between segments of 
loop. The key point of XY duality is that we can equally 
well view these loops as worldlines of charged particles 
[39] . In order to be able to construct geometrical cor- 
relation functions, we must introduce additional species 
of bosonic and fermionic particles, using a supervector 
?/? rather than a single complex scalar z [28]. A concrete 
example of this correspondence between loops and super- 



vectors was displayed in Sec. Ill A The A.J term in the 
action implies that the gauge field couples to i/? in the 
usual way. The interactions between vortex segments in 
J c translate into local interactions for 



£duai = |(V - iA)xJ>\ 2 + M |Vf + A|V| 4 + 



(37) 



The line of transitions out of the extended phase at K < 
K c and \H\ > is governed by the (compact) CP k ^ k 
model. When H = we obtain instead the NCCP k \ k 
model, a supersymmetric extension of the conventional 
dual theory for the XY transition. We will consider the 



phase diagram of this noncompact theory in Sec. IV E 
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The discussion of tricolour percolation in Sec. |III A] fo- 
cussed on the line pc = Pb in the two-dimensional pa- 
rameter space of the model. As a caricature of tricolour 
percolation away from this line, we can modify the po- 
tential of Eq. ([28]) to H cos + H 3 cos 3(9 - 0). Expand- 
ing in Hs leads to the charge three monopoles we expect 
from Sec. Ill A| (the connection between monopoles and 



anisotropy has been noted in [38 J. For general these 
have a complex fugacity i^ 3 e ±3z0 /2, just as Eq. (23 ) leads 



to complex monopole fugacities for general £>r, pq, Pb 
[82] . Changing the bare values of H3 or is not expected 
to change critical behaviour in the extended phase or at 
the CP k \ k critical point, where the low energy degrees 
of freedom are the neutral ones in the supermatrix Q. 
In particular, we expect to see the same behaviour for 
tricolour percolation everywhere on the critical line [83] . 



C. Statistics of vortex lines 

Having identified the field theory for vortex lines in short- 
range correlated w(x), we now use it to answer some 
natural questions about their statistics. 



1. Geometrical correlation functions 

The geometrical correlation functions of interest to us, 
the probabilities of various geometrically-defined events, 
can be expressed as correlation functions of i/?s. At a for- 
mal level the relationships follow from (17), as indicated 



in Section [III A 2[ but they can easily be anticipated by 
thinking of vortices as particle worldlines. 

Let G(x — y) be the probability that a single vortex 
visits both x and y. (Strictly, since we are working in the 
continuum and taking the vortices to be of infinitesimal 
thickness, we should define a small e and say 'visits re- 
gions of size e around x and y\) At the critical point and 
in the localized phase all vortices are finite, and the con- 
tributions to G(x — y) come from configurations in which 
a finite loop visits the two points. In the extended phase 
there are also contributions of a different type, with x 
and y lying on a single infinite strand. We first discuss 
the cases where these infinite strands are absent. 

The operators needed to construct geometrical corre- 
lators are combinations of the gauge-invariant bilinears 
Qap _ (^^s a P. In fact we will need only Q a ^ 

for a,/3 < k + 1 and ft/ft i.e. z a z* 13 for a ^ f3. As 
discussed in Sec. [TTJ such an insertion restricts the loop 
expansion of the partition function to configurations in 
which the point x has an incoming strand of colour a and 
an outgoing strand of colour f3 - or equivalently, configu- 
rations in which there is a loop passing through x which 
changes colour from a to /3 at x. Thus we see that the 
correlation function (Q a ^ (x)Q f3a (y)) forces a single loop 
to visit both x and y: 



Similarly, the probability that m loops connect the vicin- 
ity of x to the vicinity of y scales like (Q 12 (x) m Q 21 (y) rn ). 
As a further example, the probability that the three sites 
x, y and z all lie on the same vortex can be written 
F(x,y,z) k (Q 12 (x)Q™(y)Q 31 (z)). 



2. Vortices in the critical region 

At the critical point, the two-point function decays as 
G(x) ~ l/|x| 1+r? , where 77 is the anomalous dimension 
of Q. Moving slightly away from criticality, a finite cor- 
relation length appears, scaling as £ ~ \h — h c \~ v if h 
parametrizes the bias: in the localized (disordered) phase 
£ is a measure of the typical size of a large loop (e.g. G(x) 
decays with the Ornstein-Zernicke form \x\~ 1 e~^^) J 
and in the extended (ordered) phase, anticipating the 
next subsection, it sets the scale beyond which correla- 
tions are Brownian. In this phase Q aquires a nonzero 
expectation value (Q) ~ (h c — h)P, which gives the proba- 
bility that a given point lies within a fixed small distance 
of an infinite vortex, or that a given link in tricolour per- 
colation lies on an infinite tricord. 

A scaling argument [601 ED] relates 77 to the fractal 
dimension df of the critical vortices, and to the exponent 
r governing the number density n(l) of vortex loops of 
length I. With n(l) ~ l~ r , 



d 



f 



11-7/ 

5 — 77 



(39) 



Hyperscaling in 3D giving f3 = v(l + rj)/2. 

Precise numerical values are available from simulations 
of a lattice loop model [24] which is also in the CP k ^ k 
universality class [33] : 

d f = 2.534(9), r = 2.184(3), v = 0.9985(15). (40) 

This yields 77 = -0.068(18). 

Clearly this critical point is not accurately described by 
mean field theory, which would require v = 1/2. However 
it is interesting to consider what the mean field prediction 
for df is: surprisingly, it is not the 'trivial' value d^ F = 2, 
but the fairly accurate value d^ F = 5/2. The reason one 
might have expected d^ F = 2 is that mean field theory 
corresponds to a free field description of the critical point, 
and free field theory is usually associated with Brownian 
walks (as in the extended phase, where free field theory 
for the Goldstone modes leads to Brownian behaviour). 
Here, however, that is not the case, because the appro- 
priate free theory is formulated in terms not of if) but of 
the composite field Q: in a replica approach, 



£ free = tr(VQ) 2 



(41) 



G(x-y)<x(Q 12 (x)Q 21 (y)). 



(38) 



This leads to 77 = and a fractal dimension d^ ¥ = 5/2. 
Since 77 is small, as is typical for critical points in 3D, 
this approximation is quite accurate. The fact that df 
is close to 5/2 rather than to two is a clear symptom 
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modes (for a, /3 = 2, k + 1), 




FIG. 7. (Color online) Dominant contributions to the 
correlation functions (Q 23 (x)Q 32 (y)) ~ |x — ?/| _2 and 
(Q 12 (x)Q 21 (y)) ~ \x — y\~ x in the extended phase as de- 
scribed in the text. The orientations of the loops are not 
indicated. 



of confinement and the fact that vortices are worldlines 
not of Q, but of its 'square root', i/?. (In the context 
of quantum magnets described by CP 1 /NCCP 1 models 
[55] , a large value of n has been used as a diagnostic 
for deconfinement, i.e. the non-compactness of the U(l) 
gauge field.) The above assumes a soft spin formulation 
for Q, which we discuss further in Sec. [IV 



3. Extended phase: Goldstone modes and Brownian loops 

In three dimensions, a random walker need never return 
to his starting point, but can escape to infinity. The same 
is true of vortices in the extended phase. Of course, in a 
finite system there is no such thing as escape to infinity: 
what takes its place depends on the boundary conditions. 
The simplest and most natural choice is to allow vortices 
to terminate on the boundary. The possibility of escap- 
ing to infinity then means that a vortex strand passing 
through the origin has a finite probability of reaching the 
boundary, even in the thermodynamic limit. If the linear 
size of the system is L, then these system-spanning lines 
have an arc length of order L 2 , since their fractal dimen- 
sion is two. On the other hand, if we were to impose 
periodic boundary conditions, these lines would join up 
into loops with typical length L 3 [25j [33j [62] . 

In order to discuss geometric correlation functions in 
the extended phase we must separate out the 'infinite' 
strands. We take a finite system in which vortex lines 
can end on the boundary, and require that any line which 
does so has a specific (bosonic) colour, say a = 1. With 
this protocol the fugacity per vortex is still one, regard- 
less of whether the vortex forms a closed loop or ends on 
the boundary. In the CP k \ k language, these boundary 
conditions correspond to applying a symmetry-breaking 
field for the CP k ^ k order parameter on the boundary 
[84] . The extended phase is the ordered phase of the 
field theory, and this boundary condition fixes the di- 
rection of the order parameter to be ij) oc (1,0, ...,0). 
We may fix the gauge freedom by taking the first com- 
ponent real and positive. Fluctuations involve k (com- 
plex) bosonic Goldstone modes, 2 ,...,0 fc+1 , as well as k 
fermionic ones, which are just the original fermions x : 
if) rsj (y/l — \cj)\ 2 — |x| 2 ,0,x). To leading order in these 



Q 



11 



Q 



al la. 



Q 



Since (cj) a (x)(j)*P (y)) oc 5 a ^\x — y\ 1 we therefore have 



\x-y\ 2 



x-y\ 



The first of these gives the probability that x and y lie on 
the same finite loop, since strands of colours 2 and 3 are 
prohibited from escaping to infinity. On the other hand, 
the second of them is dominated by configurations in 
which a single infinite vortex passes through both points 
- that is, they are connected by a single finite strand, 
rather than by two strands as for points sharing a finite 
vortex. 

In general, one finds that the probability Gm(%) that 
two regions separated by a distance \x\ are connected by 
M (finite) strands of vortex scales like |x| _M : 



G M (x) ~ \x\ 



-M 



(The correlator G discussed in the context of the localized 
phase and critical point corresponds to G2; in this regime 
Gm — for odd M.) These simple exponents indicate 
that the vortices behave on large scales like free random 
walks - each factor of \j\x\ is just the probability that 
a random walker, starting in one of the regions, happens 
to visit the other. It follows that the fractal dimension is 
two, and the number density of loops of length I behaves 
as n(l) ~ r 5 / 2 . 

As another example, we easily check that the probabil- 
ity that three distant sites lie on the same vortex behaves 
as 



F(x 1 ,x 2 ,x s ) 



ri2 + r 23 + r 31 

^12^23^31 



Tnq \Xi 



u 3 h 



being dominated by contributions from infinite vortices, 
while the probability that they lie on the same finite vor- 
tex scales as (ri2r23r3i) _1 . 

This result - Brownian statistics for vortices in the 
extended phase - is expected on the basis of numerical 
simulations [3 [8j [10]. Surprisingly, a clear theoretical 
derivation has not previously been provided, although a 
heuristic explanation [8 has been given by analogy with 
polymers in the melt [26 . 



IV. VARIATIONS AND PERTURBATIONS 

So far, we have discussed universal behaviour for vor- 
tices in short-range correlated environments. In this sec- 
tion we consider the stability of this universality class 
to perturbations, including long-range correlations in the 
complex field and a perturbation leading to a crossover 
to conventional percolation, and also the universal be- 
haviour of other kinds of vortex problem: unoriented vor- 
tices in nematics, and vortices in the vicinity of the XY 
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critical point. We also briefly discuss the 6 — e expansion 
for vortices. 

A useful language for considering perturbations of the 
CP k \ k model is the soft spin Lagrangian for the gauge- 
invariant composite field Q, Eq. Since some pertur- 
bations cannot be written in the supersymmetric theory, 
and also for simplicity, we use the replica formulation in 
which Q is an n x n traceless Hermitian matrix with the 
implied limit n —> 1. 

Microscopically, this matrix is given by Q = zz^ — 1/n 
with z^z = 1, so satisfies the nonlinear constraint (Q + 
1/n) 2 = (Q + l/n) in addition to the linear one trQ = 0. 
In a soft spin approach, we might imagine imposing a 
softened version of the former via a potential 



tr((Q + l/n) 2 -(Q + l/n))' 



= const. 



n 2 — 6n + 6 



trQ 2 



2(2 -n) 



trQ 3 + tr( 



Of course these 'bare' coefficient values have no signifi- 
cance, except to draw our attention to the fact that (ex- 
cept in the case n = 2, where there is an additional sym- 
metry Q — )> — Q) there is a symmetry- allowed cubic term 
[41] [43] [63] . For n > 2, this term implies that the mean 
field prediction (expected to be valid in a large number of 
spatial dimensions) is a first order transition. However, 
this conclusion does not apply to the replica limit n — >• 1. 
Such a situation is familiar from percolation [64] . 



A. Upper critical dimension and 6 — e expansion 

The cubic term implies that the upper critical dimension 
for CP n_1 in the limit n — » 1 (and CP k \ k ) is expected 
to be six. Certain correlation functions for Anderson 
localization in symmetry class C are also described by 
this field theory, and can in fact be mapped to correlation 
functions in classical loop models described by CP k ^ k 
[24j|29> 33 > 44-46 , so we expect six to be an upper critical 
dimension for that problem too [43] . 

We may also consider the 6 — e expansion, which was 
performed for theories with a single cubic invariant in 
[4QU42] . Eq. ^ at n > 1 is precisely the theory con- 
sidered in [41], in the form given by writing Q in terms 
of (n 2 — 1) real scalar fields t{ and the SU(n) genera- 
tors Ti via Q = t{Ti. With the tensor dijk defined by 

£soft spin = 2^^) 2 + + -^dijkUtjtk. (42) 

Coefficients in the e expansion are given by contractions 
of the ds [40]. In the limit n —> 1, to order e 2 



1 5e 733e 2 



76 27436 



e 166e 2 
19 6859 



(43) 



In three dimensions, this gives v ~ 0.94 and 77 c± —0.38, 



are compatible with the 0(e) results for class C local- 
ization in [33] . Note that the sigma model formulation 
Eq. ([9| leads to an alternative approximation scheme for 
exponents via the 2 + e expansion, which has been per- 
formed to high orders - see [44] and references therein. 

The field theory in higher dimensions can again be re- 
lated to loops, and it is natural to conjecture that in d 
dimensions these have the statistics of zero lines in a ran- 
dom (d— l)-component vector field, {w\, w^-i). (Such 
lines are oriented in any number of dimensions, with ori- 
entation given by the vector e fI1 ... fJjd d fJj2 w 1 ...d fId w d - 1 .) 



B. Vortex gluing and percolation 

In cases where intersections of vortices cannot be ne- 
glected (such as in the lattice XY model), a conven- 
tion is required for defining geometrical correlation func- 
tions, or in other words for determining which vortex 
strands are regarded as connected. We may either adopt 
a convention in which vortices remain topologically one- 
dimensional, or we may 'glue' intersecting strands to- 
gether to form branching, net-like clusters [T7J [18] , and 
our choice may affect the observed critical behaviour. 
We now address this issue. The punchline is that with 
the former convention (when vortices are topologically 
one-dimensional), the results of the supersymmetric dual 
theory continue to apply even in the presence of inter- 
sections, whereas vortex gluing induces a crossover to 
conventional percolation. This crossover can be seen in 
the soft spin language. Another reason for clarifying the 
relation to percolation is that it has previously been be- 
lieved that some or all of the exponents for tricolour per- 
colation coincide with those for conventional percolation 
[6j [10] [65] . The soft spin formulation makes clear that 
this is not the case. 

A microscopic model with intersections can be got by 
modifying the lattice for tricolour percolation (specifi- 
cally, we can replace the bcc with the fee lattice, so the 
cells are rhombic dodecahedra). Alternately, we can use 
the lattice XY model. Considering such models shows 
that the dual SUSY theory {CP k \ k or NCCP k \ k ) still 
applies if we resolve vortex intersections according to 
the 'stochastic rule' used in some numerical simulations 
[TT] [T8] : at an intersection with q incoming and q outgo- 
ing vortex strands, we choose randomly from the q\ ways 
of pairing the former with the latter (Fig. [8| . This rule 
is a simple consequence of formulas like Eq. (17) for the 
traces at a single site. 



To induce a crossover to percolation, we modify the 
rule by gluing intersecting vortices with a probability 
c > [T7] [18] [85] . In the replica language, each strand 
acquires a colour a = l,...n, and to glue strands we 
must force their colours to coincide; the probability G per c 
of two links lying in the same cluster will then be given 

by 



to be compared with Eq. (40). The expressions Eq. (43) 



(Q aa Q aa ) 



(44) 
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being simply related to the probabilities of their being the 
same or different colours. We can check in the graphical 
expansion of a lattice model with intersections (we omit 
details) that the perturbation induced by the gluing is of 
the form SC ~ —c ^ a \z a \ 4 . As expected, this operator 
creates a meeting of four worldline strands of the same 
colour. The perturbation breaks rotational symmetry in 
replica space, yielding an additional mass term for the 
off-diagonal elements of Q in £ so f t spin . Subtracting a 
multiple of trQ 2 , which in the absence of the symmetry- 
breaking term would merely shift the critical point, 



SC 



\Q 



a/3,2 



(45) 



SC is expected, on the basis of its naive scaling dimen- 
sion, to be relevant at the CP k ^ k critical point. (Both SC 
and the perturbation discussed in Sec. 4C can be written 
in terms of the traceless symmetric tensor Sq] built out 
of Q a ^Q l8 , modulo a trQ 2 term.) Symmetry considera- 
tions then lead us to anticipate a crossover to behaviour 
in the universality class of the n-state Potts model in 
the limit n — >• 1, which is a well-known description of 
percolation [58]. Note that the Goldstone modes in the 
extended phase acquire a mass, consistent with the fact 
that in conventional percolation the percolating phase is 
massive. 

The above picture implies that the bias-induced vor- 
tex transition is in a distinct universality class from that 
of percolation. Previously this has been obscured by the 
fact that some exponents are similar [6][TOJ[65]. In partic- 
ular, the fractal dimension of critical tricords was found 
[6] to be df = 2.54(1) (see also Eq. (|40|), close to that 
of critical percolation clusters, for which a more recent 
study [66] gives df rc = 2.5226(1). In the light of the 
above, this similarity is not surprising, and just reflects 
the small anomalous dimension of Q at each of respec- 
tive critical points. (The result dj erc ~ 5/2 can also be 
seen without invoking field theory, via a Flory-like ar- 
gument due to de Gennes [67.) The correlation length 
exponents for the two transitions are clearly different - 
compare ^ perc = 0.8733(5) [66] with Eq. (gob. 




FIG. 8. Left: vortex strands intersecting at a site of the 
dual lattice (cube of the original lattice) in a cubic lattice XY 
model. According to the 'stochastic' rule, we choose randomly 
from the ways of resolving the intersection (centre, right). 




FIG. 9. (Color online) 'Rewiring' of vortices induced by per- 
turbing the oriented vortex problem. 



C. Unoriented vortices 

So far we have discussed oriented vortices in random com- 
plex fields, which are associated with the fundamental 
group 7Ti(6' 1 ) = Z. It is natural to expect unoriented 
vortices, which appear when the relevant fundamental 
group is equal to Z2, to correspond to a different field 
theory, since the distinction between positively and neg- 
atively charged fields (?/> and 1/?*) in the CP k ^ k model 
is due to the distinction between vortices and antivor- 

tices. To obtain the field theory for Z2 vortices from the 
cp k\k model 

we must introduce a perturbation which 
reduces the symmetry (from superunitary to orthosym- 
plectic); the introduction of crossings into 2D loop mod- 
els for non-crossing loops leads to the same symmetry 
breaking [70] . 

Unoriented vortices occur (for example) in nematic or- 
der parameters, where the order parameter manifold is 
RpN-i w ith N > 2. This is the space of TV-component 
unit vectors v with the identification v ~ — v. Vortices 
of this kind, for N = 3, have been observed in nematic 
liquid crystals [68j [69] and their fractal geometry has 
been studied numerically in the context of cosmic strings 

nana. 

As for oriented vortices, we may consider the phase di- 
agram of vortices in short-range correlated RP N ~ 1 fields 
as a function of bias favouring some point in the order 
parameter manifold. At small bias, there is an extended 
phase similar to that for oriented vortices [TOJ [13] . We 
argue below that this phase is again Brownian (thus we 
disagree with earlier suggestions that the fractal dimen- 
sion is less than two in this phase [TO]). However, this 
Brownian behaviour results from the Goldstone phase of 
a field theory different to the CP k \ k model, so the bias- 
induced critical point is in a different universality class to 
that for oriented vortices. The initial investigation of this 
transition by Strobl and Hindmarsh [TO] is probably not 
sufficiently precise to confirm this, and further numerical 
work would be interesting. 

To determine the field theory, we access the unori- 
ented vortex problem by applying a perturbation to that 
of oriented vortices. Take an RP 2 order parameter 
v — (vi,V2, vs) with v 2 = 1 and v ~ —v. Now, the 
subspace ^3 = is RP 1 = S 1 ; thus in the limit where 
fluctuations in vs are completely suppressed, we return 
to the oriented vortex problem. Turning fluctuations in 
V3 back on allows the orientation of a vortex to fluctuate 
along its length. Specifically, the leading effect of reintro- 
ducing fluctuations in vs is to allow configurations of the 
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type shown in Fig. [9j where two vortices passing close to 
each other are rewired in a manner 'incompatible' with 
the orientations of the strands. 

This perturbation is most easily discussed in the 
replica theory. In terms of wordlines, the right hand side 
of Fig. [9] corresponds to the annihilation of two particles 
of the same colour /3, and the creation of two particles 
of colour a, so the perturbation leads to the following 
SU (n)-symmetry-breaking term: 



a, (3=1 



){ z P z P) = -\z T z\ 2 . 



(46) 



In the soft-spin theory for Q, which we decompose into 
real and imaginary parts, Q = Qr + iQi, this yields 



(the subscripts give the dimensions of the unit matri- 
ces), similar to the 2D supersymmetrized O(n) model 
described in [70 in the context of lattice loop models 
with unoriented loops, but with the additional identifi- 
cation ~ — <£. Note that we can also cast the perturba- 



tion (46) in supersymmetric language. We must take the 
CP fe l "model with k = 21. Then (46) should be replaced 
with 

6C = -(^ T A^)(^ f A^*). 

It can be checked in a lattice model (such as_the mod- 
ification of tricolour percolation described in IV B) that 



this is the desired perturbation and that it gives all loops 
the correct fugacity of one. 



(47) 



a,/3=l 



We see that allowing fluctuations of the vortex orienta- 
tion yields a mass term for the imaginary part of Q. The 
simplest assumption is that this leads to a crossover to 
a similar theory but with real rather than complex Q. 
The resulting traceless symmetric matrix is the order pa- 
rameter for RP 71 ' 1 as opposed to CP 71 ' 1 . (This RP 71 ' 1 
field should not be confused with the RP N ~ X field which 
hosts the vortices.) This can be easily checked in the 
case n = 2, where the CP 71-1 model is the 0(3) model, 
and where the perturbation is a mass term for one of the 
components of the 0(3) spin, leading to a crossover to 
the XY, or equivalently RP 1 , universality class. 

So we expect that unoriented vortices in a (short-range 
correlated) random nematic field are described by the 
replica limit of the RP 71-1 sigma model. Although for 
brevity we have given only a perturbative argument, ad- 
ditional evidence comes from the existence of lattice loop 
models with RP n ~ x degrees of freedom [33]. It is also 
possible to construct a tricolour-percolation-like model 
which can be formulated as a lattice gauge theory with 
the same symmetry. 

Note that in the extended phase it makes no differ- 
ence whether we think of the field theory for unoriented 
vortices as the replica limit of an RP 71-1 sigma model, 
with target space RP 71-1 = 5' n_1 /Z2, or an 0(n) sigma 
model with target space S™ -1 . The topological difference 
between these manifolds is not seen by the Goldstone 
modes. However the difference does become important 
at the critical point: the 0(1) sigma model has a thermo- 
dynamic (Ising) transition, whereas the vortex problem 
has a geometric transition which is invisible in thermo- 
dynamic quantities. This is reflected in the fact that the 
manifold RP n ~ x becomes a point when n = 1. 

In SUSY language, the replica limit of the RP n ~ x 
model becomes a sigma model for a 'real' supervector 
field $ = (S 1 , S 2l+1 , 77 1 , rj 21 ), with <M$ = 1, where 



A = 



l 2 i+i 
1, 
-1/ 



D. Stability to long range correlations in w(x) 

We have discussed short-range correlated w(x) at length. 
How short-ranged do the correlations have to be, in order 
not to alter the universal behaviour in the extended phase 
or at the critical point? 

Weak correlations in w(x) can be addressed in a stan- 
dard way using the extended Harris criterion [58j I7T] . 
Correlations in \w\ are more relevant than correlations 
in argiu, and we will consider only the former. Letting 
these correlations decay as l/\x — y\ A at large distances, 
one finds that the extended phase is extremely robust, 
being unaffected by sufficiently weak correlations for any 
A > 0. (The statistics of vortices in correlated complex 
fields have been studied numerically in [14], and the re- 
sults appear compatible with this claim.) The critical 
theory on the other hand is stable so long as 



A > 2/v ~ 2.003(3). 



(48) 



These criteria may be derived by viewing correlations in 
w(x) as due to quenched, correlated disorder in the pa- 
rameters of the probability distribution for w [71] (for ex- 
ample, quenched disorder in the parameters £>r, pq and 
Pb of tricolour percolation) and averaging over this dis- 
order to give an effective Lagrangian with non-local cou- 
plings between the operators of the CP k ^ k model, whose 
relevance or irrelevance can then be determined. This 
average does not require use of the replica trick, since 
the partition function takes a trivial value which is inde- 
pendent of the disorder realization. 



In addition, correlations decaying slower than Eq. (48) 



demands may still count as short-range if they oscillate 
with distance [72]. This phenomenon occurs for random 
superpositions of plane waves exp (ik.x) with fixed |fc|. 
According to Berry's random wave model [73], such su- 
perpositions 



w(x) 



k 2 =E 



k a(k)e 



\k.x 



(49) 



with Gaussian, delta-correlated a(k\ give a good descrip- 
tion of the statistics of eigenfunctions of the Laplacian, 

- V 2 w = Ew, (50) 
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in chaotic d- dimensional domains. In 2D, the version 
of this problem with real w has been studied exten- 
sively in the context of quantum chaos [74 j, and the 
nodal lines w(x) = have been found to have the statis- 
tics of percolation cluster boundaries. Bogomolny and 
Schmit [72] reconciled this with the square root decay of 
the correlation function at large distances, (w(x)w(y)} ~ 
cos(E 1 ^ 2 \x — y\ — 7r/4)/^/\x — y\ - which would usually 
be slow enough to make percolation results inapplicable 
[7T] - by taking oscillations into account. A simple gen- 
eralization of this analysis shows that in 3D, for complex 
w, the collapse transition of the lines w(x) = h induced 
by increasing \h\ is expected to be in the CP k ^ k uni- 
versality class, despite the slow decay of the correlator 
(w*(x)w(y)) oc sin( J E 1 / 2 |x - y\)/\x - y\. 



E. Vortices near the XY critical point 

It has long been recognized that vortices play an impor- 
tant role in the phase transition of the 3D XY model. 
XY duality makes this more concrete, giving a dual 
formulation in which the transition is represented by 
the condensation of a 'vortex' field. A natural ques- 
tion is thus whether the thermodynamic phase transi- 
tion coincides with the geometric phase transition for 
the vortices, or whether the two transitions are sepa- 
rate. Simulations of lattice XY models [TTJ [18j [76] sug- 
gest that the geometric transition occurs slightly inside 
the XY ordered phase, but very close to the the thermo- 
dynamic transition (letting K be the coupling for an XY 
model on the cubic lattice, the most accurate determi- 
nation of the geometrical critical point in Ref. [17] gives 

(#geom _ ^therm^therm _ 4 x 1Q -4) Vortices must of 

course be defined using the stochastic rule (Sec. IV B). 

At first sight, the NCCP k \ k formulation of this prob- 
lem 



-NCCP k \ k 



|(V - iA)^\ 2 + «(V x A) 2 + /i|V>| 2 + A|^f 



implies that the transitions should coincide, occurring 
where the field i/? condenses. However, this assumption 
is not justified, because in the theory with a noncompact 
gauge field we must distinguish between the condensation 
of ?/? and the condensation of the bilinear Q. 

The geometric transition into the extended phase is sig- 
naled by the appearance of an expectation value of Q, i.e. 
by the spontaneous breaking of supersymmetry: (Q) ^ 0. 
On the other hand the thermodynamic transition into the 
XY disordered phase is signaled by the onset of the Higgs 
mechanism, due to the condensation of i/?, and the gen- 
eration of a mass for the gauge field. Following standard 
convention we denote this ^ (though of course this 
expectation value is zero by gauge invariance). In the 
simplest scenario, these events would occur simultane- 
ously. However the numerical results [17] [18] [86] lead us 
to consider the possibility of a phase with ( < 0) = but 
(Q) 0. Heuristically, this is a pair condensate [75] of 



vortices in which the XY order is not disrupted because 
vortices and antivortices are locally paired. 

If the geometric transition does occur within the XY 
ordered phase in (for example) lattice XY models, a nat- 
ural guess is that it is in the same universality class as 
the bias-induced transition. Note that the nonzero value 
of the XY order parameter yields a bias. In this sce- 
nario, long distance behaviour precisely at the geometri- 
cal phase transition would be described by a Lagrangian 
of the form 



C = «(V X A) 2 + £ cri tical CP*\4Q\, 



(51) 



where the massless noncompact gauge field A (describing 
Goldstone fluctuations in the XY order) is decoupled, at 
long wavelengths, from the neutral degrees of freedom in 
Q. These are then described by a critical compact CP k \ k 
model. The symmetry-allowed couplings between A and 
Q are irrelevant at this critical point. This scenario is 
also compatible with the extended Harris criterion (note 
that in the Goldstone phase correlations in the modulus 
of the complex field decay much faster than correlations 
in its argument). 

Previous attempts have been made to relate the frac- 
tal dimension of vortices in critical XY or Abelian Higgs 
models to local correlators in field theory using the dual- 
ity between these two theories (2TJ [76j [77] . However, the 
authors did not make use of SUSY or replica, which are 
necessary in order to form nontrivial geometrical corre- 
lators. As a result the scaling relations put forward were 
not correct, as has been previously noted [22j[78]. Note 
that vortices in Abelian Higgs models, with Lagrangian 
£ sc = |(V-m)$| 2 + £(Vxa) 2 + V(<I>), can be treated in 
a similar manner to vortices in ungauged fields, leading 
to a dual SUSY theory without a gauge field (Eq.|2| [87]. 



LOOPS IN TWO DIMENSIONS 




FIG. 10. (Color online) Percolation cluster boundaries. We 
will view the black clusters as surfaces in the graphical ex- 
pansion of a lattice gauge theory. 

We now turn to a different topic, namely 2D percola- 
tion. It is interesting to see how the relation to a lattice 
gauge theory is modified for two dimensional determin- 
istic walks in a random environment in the universality 
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class of percolation hulls (cluster boundaries). Essen- 
tially the idea will be to take the black clusters in per- 
colation (Fig. 10) to be sheets of random surface in the 
graphical expansion of a lattice gauge theory. This can 
be compared with the treatment of tricolour percolation 
in Sec. Ill A[ where the random surfaces generated by a 
lattice gauge theory were interpreted as domain walls. 
The end result in 2D is again a 

C pk\k mode ^ but this 

time with a topological term. This is a known contin- 
uum representation of percolation [29 , 32 , derived previ- 
ously by mapping bond percolation to a supersymmetric 
spin chain. The present approach makes the appearance 
of the term transparent. There is also closely related 
work mapping (non-SUSY) lattice CP 71-1 models with a 
term to solid-on-solid models on the honeycomb lattice 
[47] . That work considers a graphical expansion in terms 
of loops very similar to the one described below. How- 
ever, since it is not concerned with percolation, it uses a 
more conventional Boltzmann weight for the gauge field 
which does not yield the simple expansion in terms of 
percolation configurations that we desire. 



1. Lattice gauge theory for 2D loops 



cluster boundaries. Fixing V, we then expand the sec- 
ond product (over links). Only one term survives after 
tracing over J7, namely that which cancels the Us on the 
on the cluster boundaries, and this term comes with if) 
'worldlines' which encircle the black clusters anticlock- 
wise. The trace over if) gives the loops a fugacity: 



z = J2p b (i-p) 



W ^no. loops^total loop length 



Here B and W are the numbers of black and white 
hexagons in the configuration V. Note that our bound- 
ary conditions mean that the region outside the lattice 
is considered to be white for the purposes of drawing the 
cluster boundaries. Percolation is the case n = 1, and 
by varying n at p = 1 /2 we can consider the well-known 
0(n) loop models for various n [47] 157 ^ 158] . 

For the second approach, we expand out only the sec- 
ond product in Z, and represent the terms by drawing in 
oriented links (ji) of the lattice for each Uijip\ipj we take. 
Tracing over i/?, only configurations of mutually avoiding 
closed oriented loops survive (see also [47]). Then for 
each loop configuration C we trace over U to get a Wil- 
son loop: 



The heuristic argument of Sec. [H] can be adapted to the 
2D case, but it is just as easy to give a concrete lattice 
construction. For definiteness, let's consider site percola- 
tion on the triangular lattice. The unit cells of the trian- 
gular lattice are faces of the honeycomb lattice, which are 
coloured black or white with probabilities p and 1 — p; the 
critical point is p — 1/2. We orient the cluster bound- 
aries (hulls) so as to have the black hexagons to the left 
- see Fig. [To] 



As for tricolour percolation (Sec. Ill A), we introduce 
supervectors f ijj i = (zj, zf +/c , xj •> Xi) on the- sites i 
of the honeycomb lattice (now with arbitrary numbers 
of bosonic and fermionic components and rj^ij) = n) and 
gauge degrees of freedom Uij on the links. We also intro- 
duce a fugacity x per unit length of loop, with x = 1 for 
uncorrelated percolation. Letting F denote a hexagon, 
the partition function we need is 



z = Tr Ii f(!-p)+p II u *i\ U[ 1+xU ^l^. 



F 



CC 



(ij)eF J (ij) 



In the products over the links in a hexagon, the links 
are oriented anticlockwise. To be definite, we consider a 
finite lattice (made up of complete hexagons) with the 
topology of the disc. 

Now consider two different approaches to the graphical 
expansion of this partition function. In the first, we be- 
gin by expanding out the first product, representing the 
terms via the rule that a hexagon is coloured white if we 
take the first term and black if we take the second. Each 
term then corresponds to a percolation configuration V , 
and comes with a product of /7s along the links of the 



Zoc^n^W(C), W(C) = ( Yl U 



(52) 



\ links eC 



The Wilson loop W(C), which is evaluated using the 
Boltzmann weight Hf((^ — P) + P Yl(ij)eF Uij) for the 
gauge field, has simple properties. It is zero if C is not 
equivalent to a percolation configuration: loops must al- 
ternate in their orientations (along any line) so as to 
allow the hexagons to be consistently coloured. Then, 
W(C) oc (1 — p) w p B . There is exponential suppres- 
sion of black areas (regions encircled anticlockwise) when 
p < 1/2, and of white areas (regions encircled clockwise) 
when p > 1/2. When p = 1/2, this exponential suppres- 
sion disappears. 

These properties suggest that the continuum action 
for the gauge field A will have a topological term 
i{6 /2tt) J d 2 xE, where E = t ilv d ll A v ^ and that 6 = tt 
corresponds to p = 1/2. Particles coupled to a gauge 
field at this value of are known as 'half- asymptotic': 
in a (1 + 1)D quantum language, they can be sepa- 
rated at no energy cost, so long as they alternate (par- 
ticle/ant iparticle) in space [47|[79]. An illuminating dis- 
cussion of half-asymptotic particles is given by Shankar 
and Murthy in [79], who also mention a worldline inter- 
pretation pertinent to our considerations here. 

This topological action in fact arises from the naive 
continuum limit of the Boltzmann weight for our lattice 
gauge theory. The product of the Us around a hexagon 
gives the flux through that hexagon, 



n 

(ij)eF 



Ui, 



exp 



d 2 xE 



(53) 
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Using this, a formal expansion in E turns the part of the 
Boltzmann weight dependent only on E into 



exp 



i(6/2ir)E-p(l-p)E 2 /2 



(54) 



with = 2pn. Since E is not small, we should be suspi- 
cious of these values of the couplings. However, the fact 
that = 7r when p = 1/2 is robust, because it reflects an 
additional symmetry which is present at p — 1/2, where 
the loop ensemble is invariant under switching black and 
white hexagons. Neglecting boundary effects, this ex- 
change equates to reversing the orientation of each loop, 
or equivalent ly to complex conjugation of i\) and i$ (note 
the necessary convention is x** = ~X f° r the fermions 
[54]). This reversal of the charge of ij) is equivalent to 
a change in the sign of in the continuum description 
below. Thus the value of which corresponds to p = 1/2 
must be physically equivalent to —6, up to boundary ef- 
fects; this is the case for = it. 

Taking account of the symmetries, and dropping higher 
terms, the continuum action is 

C = i(6/27r)E + |(V - iA)^\ 2 + /i|^| 2 + A|^| 4 , (55) 

with (0-tt) ~ (p- 1/2). 

Let us briefly discuss the phase structure of this the- 
ory when n = 1 and p = 1/2, under the assumption that 
A > 0. In order to vary /i, we introduce an Ising cou- 
pling of strength J between the colours of the hexagons, 
setting x = e~ 2J in the partition function for the lattice 
gauge theory. For strong Ising coupling, J < J c , do- 
main walls are suppressed and and i\) is a massive field. 
Decreasing J, we pass through the Ising transition, at 
which i\) condenses; since this is a thermodynamic tran- 
sition as well as a geometrical one, it is nontrivial even 
when we remove the fermions by setting k = 0, leaving 
a single complex scalar z [47 j. (It is interesting that the 
Ising transition can be described either in terms of a sin- 
gle boson z coupled to the 6 = ir gauge field or - as we 
can deduce from Eq. (55) by using bosonic duality and 
then the bosonization rules - as a single Dirac fermion 
coupled to a = ir gauge field; the Ising transition in 
the latter theory has been discussed in [79] •) Finally, for 
weak Ising coupling, J < J c , we flow to the percolation 
critical point. In this regime domain walls proliferate and 
i\) is condensed. Since i/> is condensed, it is natural here 
to use a sigma model with a fixed length for i/^?/?, as in 



VI. SUMMARY 



of vortices in disordered systems, a topic that has been 
studied numerically in diverse contexts but has lacked a 
field-theoretic description. By mapping a lattice model 
for vortices in a random complex field to a lattice gauge 
theory, and also via a continuum treatment based on XY 
duality, we related geometrical correlation functions for 
vortices to correlation functions in the CP k \ k model, and 
explored simple consequences of this correspondence. For 
unoriented line defects, such as vortices in nematic fields, 
we argued that the appropriate field theory is the RP 2l \ 21 
model. These models play an important role in the classi- 
fication of universality classes of geometrical critical phe- 
nomena in three dimensions. 

We also argued that similar mappings can be fruitfully 
applied to 2D random curve ensembles, such as contour 
lines in a random height field, for which a lattice regu- 
larization is provided by site percolation and which are 
described by CP k \ k with a term. Future work will ap- 
ply similar ideas to other random curve problems in 2D. 

Although we have focussed on vortices, we know from 
the results of [33] that CP k \ k and RP k \ k models also ap- 
ply to a class of three-dimensional loop models which un- 
dergo thermodynamically trivial geometrical phase tran- 
sitions, and we expect these field theories to be generic 
descriptions for line defects that are topologically one- 
dimensional (and cannot branch, or terminate at a dan- 
gling end) or for deterministic walks in a random environ- 
ment such as trajectories in a Lorentz lattice gas. Note 
that while at first sight the loop models seem remote from 
vortices, configurations in one of them (that on Cardy's 
'L' lattice [35]) can be specified in terms of percolation 
configurations on two interpenetrating lattices, and thus 
in terms of two height fields (albeit on different lattices) 
taking values ±1. Heuristically, then, there is a relation 
to zero lines of a complex (i.e. two-component real) func- 
tion here too. 

It is interesting that the distinction between compact 
and noncompact CP n_1 models, much discussed recently 
in the context of deconfined criticality 5^ , appears even 
in the replica limit n — » 1. It remains to be seen 
whether this actually gives a new universality class for 
vortices, or whether the geometrical transition described 
by NCCP k \ k is inevitably separate from the thermody- 
namic (inverted XY) transition and reduces to t he uni - 
versality class of the compact CP k ^ k model (Sec. IV E). 
On the other hand, we do believe that unoriented vor- 
tices (or loops in an appropriate loop model) will show a 
localized - extended transition in a distinct universality 
class from that of oriented vortices - it remains to per- 
form detailed simulations for this universality class, and 
to confirm that distinct exponents are indeed obtained. 



Line-like topological defects are ubiquitous in three- 
dimensional systems, and continue to be of theoretical in- 
terest; for example, an exciting recent development is the 
realization that in certain quantum systems vortices can 
sustain topologically protected zero-energy states [80] . 
Here we have addressed the universal fractal geometry 
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Charge-three monopoles would play an important role 
in a transition into the localized phase, at zero bias, in- 
duced by Potts-like interactions between the colours of 



the cells (these can be included by modifying the Boltz- 
mann weight for the gauge field of Eq. 



13) 
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It is a matter of convention whether the charge one or the 
charge three monopoles have the non-positive fugacity. 
While a complex fugacity for monopoles may seem a 
drastic perturbation, note that this fugacity will natu- 
rally renormalize towards positive values. Crudely, after a 
coarse graining step we cannot distinguish a charge-three 
monopole from three charge-one monopoles, so the renor- 
malized fugacity ~ aH 3 e 3i@ + bH 3 + . . . with a, b > 0. 
For tricolour percolation this is done by modifying the 
Boltzmann weight for the boundary links, replacing 
Ufjiplipj with UfjZ^Zj. In the scaling limit this yields 
a boundary magnetic field having the effect described. 
In terms of connectivity, setting c = 1 is equivalent to 
the 'maximal' definition of |17| . 

While [171 118] appear to show two separate transitions, 
scaling collapse for the geometrical observables in [18] 
required the use of the standard XY correlation length 
exponent (A. M. J. Schakel and E. Bittner, private com- 
munication), which is surprising unless either it or the 
splitting of the transitions can be attributed to finite size 
effects. 

Here a is a dynamical field. The phase structure of the XY 
model in a fixed external gauge field has been considered 
in Z. Tesanovic, Phys. Rev. B 59, 6449 (1999). 



